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In  this  study,  the  spatial  variation  of  potential  environmental  impacts  of  bioenergy  crops  is  quantita¬ 
tively  assessed.  The  cultivation  of  sugar  beet  and  Miscanthus  for  bioethanol  production  in  the  North  of 
the  Netherlands  is  used  as  a  case  study.  The  environmental  impacts  included  are  greenhouse  gas  (GHG) 
emissions  (during  lifecycle  and  related  to  direct  land  use  change),  soil  quality,  water  quantity  and  quality, 
and  biodiversity.  Suitable  methods  are  selected  and  adapted  based  on  an  extensive  literature  review.  The 
spatial  variation  in  environmental  impacts  related  to  the  spatial  heterogeneity  of  the  physical  context  is 
assessed  using  Geographical  Information  System  (GIS).  The  case  study  shows  that  there  are  large  spa¬ 
tial  variations  in  environmental  impacts  of  the  introduction  of  bioenergy  crops.  Land  use  change  (LUC) 
to  sugar  beet  generally  causes  more  negative  environmental  impacts  than  LUC  to  Miscanthus.  LUC  to 
Miscanthus  could  have  positive  environmental  impacts  in  some  areas.  The  most  negative  environmental 
impacts  of  a  shift  towards  sugar  beet  and  Miscanthus  occur  in  the  western  wet  pasture  areas.  The  spatially 
combined  results  of  the  environmental  impacts  illustrate  that  there  are  several  trade  offs  between  envi¬ 
ronmental  impacts:  there  are  no  areas  were  no  negative  environmental  impacts  occur.  The  assessment 
demonstrates  a  framework  to  identify  areas  with  potential  negative  environmental  impacts  of  bioen¬ 
ergy  crop  production  and  areas  where  bioenergy  crop  production  have  little  negative  or  even  positive 
environmental  impacts. 
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1.  Introduction 

In  recent  years,  sustainability  issues  related  to  bioenergy  have 
increasingly  become  a  point  of  discussion  in  scientific,  public  and 
political  arenas.  At  several  levels,  initiatives  for  sustainability  cri¬ 
teria,  codes  of  conduct  and  protocols  have  been  developed  to  deal 
with  this  issue  [1-5].  The  steps  towards  sustainable  biomass  cer¬ 
tification  system  are  well  described  by  Lewandowski  and  Faaij  [6] 
and  van  Dam  et  al.  [7,8].  However,  sustainable  bioenergy  criteria 
are  currently  hardly  implemented  in  (inter-)  national  policies.  In 
addition,  initiatives  regarding  formulation  of  codes  of  conduct  are 
mainly  developed  in  a  general  and  top  down  approach  and  do  not 
account  for  the  significance  of  regional  variation  in  adoption  and 
implementation. 

In  the  scientific  arena,  sustainability  of  bioenergy  has  been  an 
important  research  objective.  In  Blottnitz  and  Curran  [9]  several 
studies  regarding  environmental  impacts  of  bioethanol  production 
chains  have  been  reviewed.  In  most  of  the  studies,  energy  balances 
and  greenhouse  gas  (GHG)  balances  of  production  chains  have  been 
assessed.  In  others,  such  as  Kaltschmitt  et  al.  [10],  Lewandowski 
and  Heinz  [11]  Hamelink  and  van  de  Broek  [12],  Malca  and  Freire 
[13],  Blottnitz  and  Curran  [9]  and  Styles  and  Jones  [14]  more  com¬ 
plete  Life  Cycle  Assessments  (LCA)  have  been  performed.  Van  den 
Broek  et  al.  [15],  Smeets  et  al.  [16,17]  Rowe  et  al.  [18]  and  Bor- 
jesson  [19,20]  are  examples  of  studies  on  the  sustainability  of 
several  bioenergy  crops  in  specific  countries  that  included  a  broad 
scope  of  environmental  impacts  and  some  even  includes  socio¬ 
economic  impacts.  The  national  scope  of  these  latter  assessments 
implies  a  relative  aggregated  level  of  analysis  without  differenti¬ 
ating  for  region  specific  conditions.  Physical  and  socio-economic 
conditions  often  vary  strongly  between  and  within  regions  and 
determine  the  design  of  bioenergy  chains:  the  chains  need  to  be 
physically  feasible  in  the  geographical  region,  but  also  need  to  be 
compatible  with  the  socio-economic  context.  In  addition,  the  inter¬ 
actions  between  the  bioenergy  chain  and  the  regional  conditions 
have  a  strong  influence  on  the  actual  environmental  and  socio¬ 
economic  performance  of  bioenergy  chains.  A  national  approach 
could  therefore  give  a  distorted  impression  of  the  sustainability  of 
bioenergy  chains  within  regions.  For  that  reason,  the  regional  level 
is  a  more  suitable  level  of  analysis  of  the  impact  of  biomass  chains 
in  order  to  safeguard  sustainability  of  bioenergy  chains  for  certi¬ 
fication  purposes.  This  is  illustrated  by  van  Dam  [21]  who  made 
a  region  specific  assessment  of  impacts  of  bioenergy  production 
in  the  Pampa  region  in  Argentina.  Although  this  study  took  dif¬ 
ferences  in  physical  context  into  account,  the  impacts  were  not 
assessed  spatially  explicit.  Regional  averages  on  physical  param¬ 
eters  were  used,  e.g.  for  soil  quality  or  water  availability.  As  the 
spatial  variation  in  determining  parameters  of  impacts  of  biomass 
chains  is  significant,  impacts  would  preferably  assessed  spatially 
explicitly. 

The  aim  of  this  paper  is  to  develop  a  methodological  framework 
to  assess  the  environmental  impacts  of  regional  biomass  chains 
in  a  spatially  explicit  and  integrated  way.  The  production  of 
ethanol  from  Miscanthus  ( Miscanthus  x  Giganteus)  and  sugar  beet 


( Beta  vulgaris  L.)  in  the  North  of  the  Netherlands  (the  province  of 
Groningen,  Friesland  and  Drenthe)  is  used  as  case  study.  The  case 
area  was  selected  because  of  the  high  data  availability.  This  case 
study  was  also  used  to  analyse  the  economic  viability  of  regional 
biomass  chains  [22].  This  study  will  explore  the  spatial  variability 
of  environmental  impacts  of  bioenergy  chains  and  builds  on  the 
previous  analysis. 

The  environmental  impacts  identified  by  several  (inter-) 
national  initiatives  regarding  sustainability  criteria  of  bioenergy 
production  were  used  as  a  starting  point  for  a  selection  of  rele¬ 
vant  impacts  and  suitable  indicators,  as  these  criteria  reflect  the 
areas  of  concern  for  large  scale  bioenergy  application  [23-25].  The 
environmental  impacts  taken  into  account  in  this  study  are:  GHG 
emissions  (both  during  lifecycle  and  due  to  land  use  changes),  soil 
quantity  and  quality,  water  use  and  water  quality,  and  biodiver¬ 
sity.  The  selection  and  adaptation  of  methodologies  to  assess  the 
environmental  impacts  is  based  on  an  extensive  literature  review 
on  impact  assessments.  In  this  study,  it  is  assumed  that  there  is  no 
competition  with  food  production  as  the  land  for  bioenergy  crops 
is  limited  to  the  arable  land  that  could  become  available  by  means 
of  more  efficient  agricultural  production.  The  developments  in  land 
availability  for  bioenergy  crops  in  the  Netherlands  is  based  on  the 
study  of  de  Wit  et  al.  [26]. 

2.  Methods 

The  level  of  spatial  detail  of  the  analysis  depends  on  the  avail¬ 
ability  and  resolution  of  the  required  input  data.  For  several 
impacts,  a  postal  code  area  was  selected  as  the  level  of  spa¬ 
tial  analysis.  For  other  impacts,  a  grid  level  of  0.1-1  km2  was 
selected.  Because  of  the  heterogeneity  of  the  criteria,  it  is  nei¬ 
ther  possible  nor  appropriate  to  draw  generic  system  boundaries 
and  select  a  single  functional  unit.  Therefore,  system  bound¬ 
aries  and  functional  units  were  decided  on  for  every  individual 
indicator. 

In  this  study,  it  is  assumed  that  only  agricultural  land  (both 
arable  land  and  pastures)  is  used  for  the  cultivation  of  bioenergy 
crops.  The  effects  of  land  use  changes  (LUC)  for  energy  crop  produc¬ 
tion  were  compared  with  the  impacts  of  current  agricultural  land 
use.  Nature  areas  were  excluded  beforehand.  This  is  in  line  with  the 
EU  directive  on  the  promotion  of  the  use  of  energy  from  renewable 
sources  [23].  The  availability  of  agricultural  land  is  therefore 
directly  related  to  the  intensification  of  agricultural  land.  As  in 
the  Netherlands  the  current  state  of  agriculture  is  very  advanced 
and  consequently  yield  gaps  are  relatively  small,  only  a  small 
proportion  of  agricultural  land  is  expected  to  become  available  for 
energy  crops  [22,26],  Since  current  fertiliser  and  pesticide  applica¬ 
tion  levels  in  the  Netherlands  already  meet  legislative  limits,  it  is 
assumed  that  the  slight  production  growth  of  the  remaining  agri¬ 
cultural  land  are  not  a  result  of  increased  inputs  but  from  improved 
management  only.  It  is  therefore  assumed  that  this  intensification 
will  not  result  in  significantly  more  environmental  impact.  Con¬ 
sequently,  the  environmental  impacts  related  to  the  productivity 
increase  of  the  remaining  agricultural  land  were  not  included  in 
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this  assessment.  The  impacts  of  cultivation  of  Miscanthus  and 
sugar  beet  were  assessed  for  the  entire  agricultural  area.  This 
allows  a  straightforward  interpretation  of  the  analysis  and  pro¬ 
vides  a  rationale  for  the  selection  of  most  appropriate  areas  energy 
crop  production. 

Only  potential  direct  effects  of  implementation  of  bioenergy 
chains  were  assessed.  Second  order  effects  and  feedback  loops,  e.g. 
the  effect  of  changes  in  nutrient  leaching  on  biodiversity,  are  hard 
to  quantify  and  to  allocate  and  are  due  to  their  complexity,  beyond 
the  scope  of  this  study.  Indirect  effects  and  displacement  effects 
due  to  LUC  were  tackled  beforehand  by  assuming  that  only  surplus 
agricultural  land  could  be  dedicated  to  bioenergy  crops.  Based  on 
the  first  order  results  of  the  assessment  of  the  spatial  variability 
of  environmental  impacts  combined  with  results  of  the  study  on 
the  spatial  variation  in  economic  viability  of  bioenergy  chains  [22], 
areas  can  be  identified  where  bioenergy  crop  production  is  eco¬ 
nomically  viable  and  have  the  least  negative  or  possibly  positive 
environmental  impacts. 

The  environmental  impact  analysis  of  LUC  to  bioenergy  crops 
on  a  spatially  explicit  and  detailed  level  requires  large  amount  of 
spatial  data  of  soil  characteristics,  climate,  current  land  use,  yield, 
crop  management,  fertiliser  and  manure  inputs  and  deposition  lev¬ 
els,  see,  e.g.  [16,21].  To  deal  with  the  large  amounts  of  spatial  data, 
a  framework  was  developed  to  be  able  to  calculate  the  impacts 
included  in  this  study. 

The  Miterra  model  used  in  our  study  is  derived  from  the 
MITERRA-Europe  model,  which  was  developed  to  assess  the  effects 
and  interactions  of  policies  and  measures  in  agriculture  on  nitrogen 
(N)  losses  on  a  regional  level  in  EU-27  [27].  Miterra  simulates  the  N 
and  phosphorous  (P)  balance,  emissions  of  NH3,  N20,  NO*,  and  CH4, 
leaching  of  N,  N03_  concentration  in  groundwater  and  changes  in 
soil  and  biomass  carbon  stocks  at  a  postal  code  level.  For  the  appli¬ 
cation  in  this  study,  the  model  was  adapted  and  input  data  and 
parameter  values  were  made  region  specific.  The  main  input  data 
of  the  model  are  crop  areas  and  livestock  numbers,  which  were 
derived  from  the  Geographical  Information  system  Agrarian  Busi¬ 
nesses  database  (GIAB)  on  postal  code  level.  For  the  assessment  of 
the  impacts  not  included  in  the  Miterra  model,  appropriate  meth¬ 
ods  were  selected  and  adapted  for  spatial  explicit  application.  In 
the  following  sections,  the  methods  used  to  asses  every  individual 
impact  are  discussed. 

2.2.  GHG  emissions 

One  of  the  main  drivers  of  bioenergy  is  the  reduction  of  GHG 
emissions.  Therefore,  the  GHG  emission  related  to  bioenergy  pro¬ 
duction  is  an  important  indicator  to  take  into  account.  The  GHG 
emission  of  the  bioenergy  chains  was  calculated  by  combining  the 
GHG  emissions  related  to  direct  LUC  and  the  GHG  emissions  during 
all  stages  of  the  ethanol  production  chain. 

2.2.2.  GHG  during  life  cycle 

GHG  is  emitted  during  the  lifecycle  of  biofuels  and  is  related  to 
cultivation  of  the  energy  crop  (diesel  for  agricultural  machinery, 
seed  pesticide  and  fertiliser  production  and  fertiliser  application), 
storage  and  transport  (heat  for  drying,  diesel  for  transport,  dry  mat¬ 
ter  loss)  and  processing  (energy  and  chemical  inputs  for  conversion 
process).  The  GHG  emissions  related  to  the  production  of  bioen¬ 
ergy  were  calculated  using  a  LCA  approach  which  is  often  applied 
on  bioenergy  production  chains.  The  ISO  14040-14049  guidelines 
articulate  a  preference  for  the  use  of  system  expansion  in  order  to 
account  for  co  products.  As  for  legislative  purposes  energetic  value 
based  allocation  is  preferred  in  the  EU  directive  on  biofuels  [28], 
this  was  applied  in  this  study.  The  calculation  of  the  GHG  emissions 
was  based  on  data  provided  byJEC  [29]  on  the  production  of  biofu¬ 
els.  For  cultivation  parameters,  universal  data  were  substituted  by 


region  specific  data  provided  in  van  der  Hilst  et  al.  [22].  In  Annex 
1.1,  the  data  used  to  calculate  the  primary  energy  requirements  and 
the  GHG  balances  are  depicted. 

2.2.2.  GHG  emissions  due  to  LUC 

GHG  emissions  due  to  LUC  are  caused  by  changes  in  soil  car¬ 
bon  stocks,  above  and  below  ground  biomass  and  residues.  When 
land  is  converted  from  one  land  use  to  another,  carbon  can  accu¬ 
mulate  (carbon  sequestration)  or  diminish  (carbon  emissions).  In 
addition,  LUC  causes  changes  in  N20  emissions  due  to  changes  in 
fertiliser  and  manure  application  and  drainage  of  organic  soils.  As  it 
is  assumed  that  there  are  no  changes  in  livestock,  livestock  related 
N20  and  methane  emissions  are  assumed  to  remain  constant  when 
land  is  converted  to  energy  crops.  The  Miterra  model  used  here 
makes  use  of  the  methodology  proposed  in  the  IPCC  guidelines  to 
calculate  the  GHG  emissions  due  to  LUC  [31].  In  Annex  1.1,  the 
parameters  used  to  calculate  changes  in  GHG  emissions  due  to  LUC 
are  presented. 

Direct  N20  emissions  from  managed  soils  related  to  different 
N  sources  (manure,  grazing,  mineral  fertiliser,  crop  residues  and 
cultivation  of  organic  soils)  and  indirect  N20  emissions  due  to  N 
leaching  and  N  deposition  were  calculated  by  Miterra  using  emis¬ 
sion  factors  of  the  IPCC  guidelines  [31,32].  The  N  surplus  was 
calculated  from  the  total  N  input  (manure,  mineral  fertiliser,  depo¬ 
sition  and  N  fixation),  the  removal  via  crop  harvesting  and  the  N 
losses  by  gaseous  emissions  and  runoff.  Part  of  the  N  surplus  is 
leached  which  results  in  direct  N20  emissions  [27].  In  Annex  1.1  and 
1.3,  the  parameters  to  calculate  the  changes  in  N20  emissions  are 
presented. 

Due  to  lower  nutrient  requirements,  less  manure  will  be  applied 
on  herbaceous  energy  crops.  As  no  change  in  manure  production 
is  assumed,  a  manure  surplus  could  be  developed.  This  potential 
manure  surplus  is  not  considered  here  as  it  is  outside  the  scope  of 
this  study.  The  annual  changes  in  SOC,  the  change  in  carbon  flux 
from  organic  soils  and  the  changes  in  annual  N20  emissions  were 
combined  to  calculate  the  changes  in  annual  GHG  emission. 

2.2.  Soil  quality 

As  soil  quality  is  a  broad  and  wide-ranging  concept,  it  was  nar¬ 
rowed  down  to  indicators  that  are  relevant  for  the  region  and  for 
measuring  the  effect  of  bioenergy  crop  production  on  the  soil.  Soil 
organic  matter  (SOM)  and  soil  erosion  hazard  are  the  indicators  of 
soil  quality  that  were  selected  in  this  study  in  line  with  the  studies 
of  Smeets  et  al.  [33]  and  van  Dam  et  al.  [21,34]. 

2.2. 2 .  Soil  organic  carbon 

As  the  SOC  content  is  closely  related  to  the  amount  of  SOM, 
SOC  is  used  as  a  proxy  indicator.  In  addition  to  the  benefits  of  car¬ 
bon  sequestration  in  the  soil  (see  Section  2.1.2),  SOC  is  also  linked 
to  other  important  functions  of  the  soil  like  water  holding  capac¬ 
ity,  nutrient  retention  and  soil  structure  [18,35].  SOC  is  therefore 
considered  to  be  the  most  prominent  indicator  for  soil  quality  [36]. 

SOC  is  affected  by  soil  type,  climate,  past  and  present  land  use, 
soil  and  water  management  [37].  In  the  Netherlands,  most  soils 
have  a  high  SOC  content.  Pastures  usually  contain  more  SOC  than 
arable  land  due  to  continuous  land  cover,  higher  root  biomass,  high 
animal  manure  application  rates  and  absence  of  tillage.  Sandy  soils 
have  generally  lower  organic  matter  content  than  clay  soils  because 
of  higher  organic  decomposition  rates  [37].  As  wet  conditions  slow 
down  decomposition  rate,  soils  with  high  water  tables  contain  gen¬ 
erally  more  SOC  [37].  Peat  lands  store  the  largest  amount  of  SOC  but 
annual  losses  due  to  cultivation  and  drainage  are  high  [37].  Changes 
in  SOC  were  evaluated  for  the  LUC  related  GHG.  The  methodology 
applied  is  discussed  in  Section  2.1.2. 
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Table  1 

parameters  of  the  WEQ.  equation. 


Parameter 

Unit 

See  Annex 

Ref. 

E 

I 

Erosion 

Soil  erodibility  index 

Potential  average  annual  soil  loss 

The  potential  annual  soil  loss  per  hectare  for  an  open  unsheltered 

ton  ha"4  y 
ton  ha-1  y 

Table  Al  .3 

a,b 

I< 

Soil  surface  roughness  factor 

isolated  levelled  field  with  a  smooth  bare  crustless  surface.  Related  to 
texture,  organic  matter,  calcium  carbonate  concentration 

Reduction  of  the  potential  wind  erosion  due  to  ridges  of  the  soil. 

factor 

Table  Al  .5 

b 

C 

Climate  factor 

Related  to  ridges  and  cloddiness  made  by  field  operations 

Local  climatic  erosivity.  Related  to  wind  speed,  precipitation  and 

factor 

Table  A1.4 

c 

L 

V 

Length  of  field 

Vegetation  factor 

temperature 

Unsheltered  distance  across  a  field  along  the  prevailing  wind  direction 
Equivalent  vegetative  cover.  Kind,  amount  and  orientation  of  the 

m 

ton  ha-1  y 

Table  Al  .5 

d 

vegetation  expressed  in  Small  Grain  equivalents  (SGe) 

a  USDA  et  al.  (2002). 
b  van  Kerckhoven  et  al.  (2009). 
c  Heijboer  et  al.  (2002). 
d  NRCS  (2006). 


2.2.2.  Soil  erosion 

In  areas  prone  to  soil  erosion,  LUC  may  have  a  significant  effect 
on  the  actual  erosion  risk.  The  main  on-site  problem  caused  by  ero¬ 
sion  on  agricultural  land  is  the  loss  of  fertile  top  soil  which  leads 
to  degradation  of  arable  soils  [37-40]  and  crop  damage  caused  by 
abrasion  or  burial  of  seedlings  or  plants  and  the  exposure  and  loss 
of  seed  [41  ].  In  addition,  the  transport  of  minerals,  organic  matter, 
residues  and  pesticides  could  cause  contamination  of  surround¬ 
ing  surface  water  [39,42].  Furthermore,  airborne  particles  due  to 
wind  erosion  could  affect  human  and  animal  health,  machinery  and 
infrastructure  [39,40].  The  main  factors  determining  the  actual  ero¬ 
sion  risk,  are  the  soil  characteristics  (especially  soil  moisture  and 
soil  structure),  vegetation  (soil  cover)  and  slope  (in  case  of  water 
erosion)  [37,38].  It  is  therefore  important  to  assess  the  effect  of  LUC 
or  farm  management  on  these  factors  in  regions  with  high  erodible 
soils  and  periods  with  high  erosive  rainfall  or  wind. 

Most  studies  regarding  erosion  focus  on  erosion  caused  by 
water,  as  this  is  the  dominant  form  of  erosion  in  most  areas  [38]. 
In  the  Netherlands,  soil  erosion  due  to  runoff  is  mainly  occurring  in 
the  sloping  areas  in  the  South-East  [37,43].  The  erosion  hazard  map 
of  the  Netherlands  by  Eppink  and  Spaan  [43]  shows  no  erosion  risk 
due  to  rainfall  runoff  in  the  North  of  the  Netherlands  as  this  area 
has  an  outspoken  flat  morphology  and  the  soils  have  a  relatively 
high  infiltration  capacity.  For  this  reason,  rainfall  related  erosion 
is  not  included  in  this  study.  Wind  erosion,  however,  does  play  a 
role  in  this  region.  Areas  underlain  by  sandy  and  reclaimed  cut-over 
peaty  soils  with  frequent  tillage  on  large  fields  without  barriers  are 
susceptible  to  wind  erosion  [37,40,43,44].  It  is  estimated  that  the 
vulnerable  area  of  wind  erosion  in  the  in  the  North  Netherlands 
amounts  up  to  60  kha. 

Several  methods  have  been  developed  to  model  wind  erosion  for 
different  temporal  and  spatial  scales,  functionalities  (specific  cir¬ 
cumstances)  and  impacts  (soil  loss,  particle  concentration).  Ideally, 
soil  erosion  is  continuously  measured  in  the  field  or  estimated  by 
very  exact  wind  erosion  models  fed  with  continues  data  on  physical 
parameters.  However,  this  is  very  time  and  capital  intensive.  In  this 
study,  the  wind  erosion  equation  (WEQ)  method  is  applied.  This  is 
a  relatively  simple  methods  that  requires  less  detailed  input  data 
and  can  be  applied  on  a  regional  level  as  demonstrated  by  Van  Ker¬ 
ckhoven  [42].  The  WEQ  estimates  the  average  soil  loss  due  to  wind 
erosion  (in  ton  ha-1  y-1 ).  Eq.  (1 )  and  Table  1 ,  give  the  parameters  of 
the  WEQ  derived  from  USDA  and  NRCS  [40],  Morgen  [45]  and  van 
Kerckhoven  [42].  The  complex  relations  between  the  parameters 
and  the  input  data  to  calculate  the  soil  loss  are  explained  in  Annex 

1.2. 

E=  /  (IKVCL)  (1) 


When  land  is  converted  from  conventional  use  to  bioenergy 
crops,  most  factors  included  in  the  WEQ  remain  constant.  The  soil 
surface  roughness  factor  could  change  slightly  (±20%)  due  to  mod¬ 
ifications  in  tillage  and  planting  practices.  The  main  factor  that  will 
change  is  the  equivalent  vegetative  cover.  The  erosion  risk  is  deter¬ 
mined  for  every  individual  crop  for  every  month  of  the  year.  To 
calculate  the  erosion  risk  for  rotations,  the  risk  of  the  rotation  crops 
are  combined  in  a  weighted  summation  based  on  the  portions  of 
the  crops  in  the  rotation.  For  grass,  it  is  assumed  that  pastures  are 
reseeded  every  10  years.  For  Miscanthus  a  lifetime  of  20  years  is 
assumed  (1  year  of  establishment,  19  years  of  full  crop  cover,  early 
spring  harvesting).  All  parameters  used  for  the  calculations  of  wind 
erosion  are  depicted  in  Annex  1.2.  In  line  with  the  WEQguidelines, 
the  most  critical  month  has  been  selected  based  on  calculations 
(see  Fig.  A2.3  and  A2.4  Annex). 

2.3.  Water  quantity  and  water  quality 

Water  use  and  water  quality  are  often  addressed  in  proposed 
sustainability  criteria  for  bioenergy  [3,24,46].  As  a  relatively  small 
proportion  of  the  total  water  consumption  for  biofuel  production  is 
used  during  the  processing  stage  (<1 0%)[47-50],  only  the  water  use 
related  to  crop  production  is  taken  into  account.  The  assessment  of 
the  impact  on  water  quality  is  limited  to  evaluation  of  the  effect  of 
change  in  use  of  agrochemicals  like  fertilisers  on  the  water  quality. 

2.3.1.  Water  use 

The  change  from  current  land  use  to  energy  crops  may  change 
the  water  balance  of  an  area  due  to  changes  in  evapotranspiration, 
runoff  and  percolation  [33  ].  The  amount  of  water  lost  through  evap¬ 
otranspiration  depends  on  crop  type,  growth  stage,  climate,  soil 
characteristics  growing  period  and  agronomic  practice  [48,50-52]. 
As  in  the  North  of  the  Netherlands  only  a  minor  proportion  of  the 
agricultural  land  is  irrigated  (1%  in  wet  years  to  11%  in  dry  years 
[53]),  water  withdrawal  for  irrigation  is  not  incorporated  in  this 
study. 

Preferably,  the  impact  of  energy  crops  on  fresh  water  availabil¬ 
ity  for  other  functions  is  assessed  on  a  water  basin  level  [50].  This 
approach  however,  requires  detailed  knowledge  and  data  about  the 
hydrologic  flows  within  a  specific  water  basin.  In  this  study,  the 
changes  in  the  water  balance  when  current  land  use  is  converted 
to  bioenergy  crops  was  assessed  on  a  spatial  level  of  1  km2  grids. 
In  addition,  a  comparison  was  made  between  the  region  specific 
water  use  efficiency  (WUE)  of  the  two  bioenergy  crops. 

To  assess  the  potential  water  depletion  due  to  the  introduction 
of  bioenergy  crops,  a  simple  water  balance  was  made  by  compar¬ 
ing  the  evapotranspiration  to  the  effective  precipitation  like  done 
in  the  studies  of  Smeets  and  Faaij  [17]  and  van  Dam  et  al.  [21]  (see 
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equation  in  Appendix  Al  .3 ).  The  calculation  method  of  water  deple¬ 
tion  is  made  time  and  spatially  explicit  by  making  use  of  spatial 
data  on  long  term  averages  of  monthly  precipitation  (1  km2)  of  the 
Royal  Netherlands  institute  for  Meteorology  (KNMI).  The  effective 
precipitation  (EP)  is  defined  as  the  rainfall  that  is  useful  or  usable 
in  any  phase  of  the  crop  production  [54]  and  is  derived  from  the 
actual  rainfall  making  use  of  the  USDA  formula  in  the  CROPWAT 
8.0  model  [55].  The  spatial  distribution  of  monthly  reference  evap- 
otranspiration  (ET0)  was  constructed  on  a  0.1  km2  grid  level  for 
the  North  of  the  Netherlands  by  means  of  aggregation  and  spa¬ 
tial  interpolation  of  daily  ET0  data  of  the  KNMI  weather  stations  in 
the  region.  Although  the  Penman-Monteith  equation  is  the  most 
widely  applied  method  to  calculate  ET0,  the  KNMI  uses  the  Makkin 
method  for  practical  reasons  [56,57]. 

The  crop  evapotranspiration  coefficients  (Kc  factors)  mainly 
depend  on  crop  type,  growth  stage  of  crop  and  climate  [58].  The 
crop  specific  decadal  Kc  factors  for  the  Dutch  situation  were  derived 
from  the  KNMI  [59].  The  spatial  0.1  km2  grid  maps  of  monthly 
ET0  were  converted  to  land  use  specific  monthly  evapotranspira¬ 
tion  maps,  using  the  crop  specific  monthly  Kc  factors.  The  average 
monthly  evapotranspiration  of  rotations  were  calculated  using  a 
weighted  summation  of  the  individual  crops. 

As  water  deficits  only  occur  during  the  summer  in  this  area  [59], 
the  assessment  of  the  change  in  total  water  deficit  focussed  on  the 
months  April-August.  The  parameters  used  for  the  calculation  of 
the  spatially  explicit  water  deficit  are  depicted  in  Annex  1.3. 

The  Water  Use  Efficiency  (WUE)  indicator  is  frequently  applied 
in  bioenergy  related  studies,  like  in  Berndes  et  al.  [48],  Dorn- 
burg  et  al.  [50],  van  Dam  et  al.  [21],  Fraiture  and  Berndes  [47] 
and  Smeets  and  Faaij  [17].  In  this  study,  the  WUE  is  used  as  a 
second  indicator  of  water  consumption.  It  provides  an  indication 
about  the  water  requirements  per  unit  crop  produced,  whereas  the 
water  deficit  methodology  only  provides  figures  for  water  use  per 
hectare.  In  order  to  assess  the  spatial  explicit  WUE  of  bioenergy 
crops,  knowledge  on  water  availability  (precipitation  and  ground 
water)  and  the  effect  of  different  growth  limitation  factors  on  evap¬ 
otranspiration  rates  should  be  considered.  As  here  are  too  many 
uncertainties  in  these  parameters,  only  a  regional  average  WUE  was 
calculated. 

2.3.2.  Water  quality 

The  use  of  agricultural  chemicals  like  pesticides  and  fertilisers 
may  contaminate  ground  and  surface  water.  A  change  in  the  use  of 
these  chemicals  due  to  LUC  in  favour  of  energy  crops  is  expected 
to  have  impact  on  the  water  quality.  In  this  study,  we  focused  on 
contamination  of  water  due  to  fertiliser  and  manure  use  only,  as 
these  cause  the  most  persistent  problems  in  The  Netherlands.  An 
oversupply  of  N  is  likely  to  cause  leaching  of  N  to  ground  and  sur¬ 
face  water.  The  N  surplus  was  calculated  in  the  Miterra  model  as 
described  in  introduction  section  of  the  methods.  Part  of  the  sur¬ 
plus  is  denitrified  and  part  will  leach  the  ground  and  surface  water. 
The  amount  of  N  leaching  was  calculated  by  multiplying  the  N  sur¬ 
plus  by  a  leaching  fraction,  which  was  based  on  the  soil  type,  land 
use  and  groundwater  level  [60].  The  N03~  concentration  in  ground 
water  was  calculated  by  dividing  the  N  leaching  flux  by  the  pre¬ 
cipitation  surplus.  For  P  there  is  no  direct  relation  between  the 
P  accumulation  in  the  soil  and  the  concentration  in  ground  and 
surface  water.  However,  when  P  application  levels  exceed  the  P 
removal  by  crop  uptake  over  longer  periods  of  time,  the  risk  of 
P  leaching  to  water  will  increase.  The  P  balance  is  calculated  by 
Miterra  making  use  of  data  on  P  inputs  (manure,  grazing,  mineral 
fertiliser  application  and  deposition)  and  P  uptake  by  crops  [32]. 
The  input  data  required  for  the  calculations  of  changes  in  N03_  con¬ 
centration  in  groundwater  and  change  in  P  balances  are  included 
in  Annex  1.3. 


2.4.  Biodiversity 

In  most  of  the  developed  sustainability  criteria  for  bioenergy, 
impacts  on  biodiversity  have  been  identified  as  an  important  area 
of  concern.  LUC  is  a  strong  driver  of  changes  in  biodiversity  [61-64]. 
Because  of  the  loss,  modification  and  fragmentation  of  habitats,  the 
(indirect)  expansion  of  agricultural  land  for  energy  crop  production 
is  perceived  to  be  a  major  threat  for  biodiversity  [65].  Biodiversity 
is  also  affected  by  LUC  related  depletion,  degradation  and  pollu¬ 
tion  of  ecosystems  and  invasive  species  [64,66].  Changes  in  habitats 
due  to  energy  crop  production  are  most  significant  when  natu¬ 
ral  areas  are  converted  to  (intensive)  agriculture  areas  [67].  Most 
bioenergy  sustainability  criteria  deal  with  this  issue  by  proposing 
process  indicators,  e.g.  by  referring  to  national  regulations  and  by 
excluding  protected  areas  and  land  identified  as  area  with  high 
biodiversity  from  bioenergy  production  [28].  In  this  study,  these 
areas  are  excluded  beforehand  and  only  changes  in  agricultural 
land  use  are  assessed.  Currently,  few  guidelines  are  available  about 
quantitative  result  indicators  and  methods  to  assess  the  impacts  of 
energy  crop  production  on  biodiversity  [46].  This  is  especially  true 
for  assessing  agro-biodiversity  on  a  regional  level. 

Over  the  last  decades,  several  indicator  systems  have  been 
developed  to  assess  changes  in  biodiversity.  These  indicator  sys¬ 
tems  vary  to  a  great  extent  according  to  scale  (global,  national, 
regional  or  local),  purpose  (policy  targets),  and  focal  area  of  biodi¬ 
versity  (species,  genetic  variation,  population  size  or  ecosystems). 
The  impact  of  energy  crop  cultivation  on  biodiversity  depends  on 
both  local  scale  effects  (choice  of  crop,  management  intensity,  veg¬ 
etation  structure,  substituted  land  use)  and  landscape  scale  effects 
(geographical  location,  scale  and  distribution  of  crops)  [68].  The 
objective  in  this  study  is  to  quantify  the  impact  of  bioenergy  crop 
production  on  agro-biodiversity  on  a  local  scale  and  identify  the 
areas  with  a  high  risk  on  biodiversity  loss.  Two  indicator  systems, 
the  High  Nature  Value  (HNV)  and  the  Mean  Species  Abundance 
(MSA)  were  selected  for  this  assessment. 

HNV  has  become  one  of  the  main  indicators  to  assess  the  impact 
on  biodiversity  of  the  2007-201 3  Rural  Development  Programs  [69] 
and  for  the  integration  of  environmental  concerns  into  the  Com¬ 
mon  Agricultural  Policy  [67].  High  Nature  Value  (HNV)  farmland 
are  certain  types  of  farmland  that,  because  of  their  characteristics, 
can  be  expected  to  support  high  levels  of  biodiversity  of  conser¬ 
vation  concern  [69].  JRC  and  EEA  developed  a  method  to  identify 
HNV  farmland  [70].  Elbersen  and  van  Eupen  [71  ]  adapted  the  HNV 
method  for  the  Netherlands  based  on  more  detailed  spatial  data 
on  HNV  parameters.  The  HNV  indicators  used  in  this  study  are 
explained  in  Annex  1.4.  Based  on  the  characteristics  of  the  HNV 
types  and  characteristics  of  the  energy  crops,  the  impact  of  a  shift 
to  bioenergy  crop  production  on  biodiversity  in  a  HNV  area  was 
qualitatively  assessed  on  a  detailed  geographical  level.  The  risk 
level  of  biodiversity  loss  was  identified  based  on  the  extent  to 
which  the  new  land  use  function  could  meet  the  habitat  functions 
requirements  of  the  important  species  identified  by  the  HNV  clus¬ 
ter  indicators  (see  Table  Al  .9).  The  areas  where  the  implementation 
of  Miscanthus  or  sugar  beet  could  cause  a  ‘high  risk  of  biodiversity 
loss’,  ‘considerable  risk  of  a  biodiversity  loss’  and  ‘both  positive  and 
negative  effects  on  biodiversity’  were  mapped.  The  HNV  approach 
does  not  provide  information  about  the  impact  of  a  shift  towards 
bioenergy  crops  on  biodiversity  outside  the  HNV  areas.  In  addition, 
the  HNV  approach  does  not  measure  the  change  in  biodiversity 
quantitatively. 

The  Mean  Species  Abundance  (MSA)  is  a  quantitative  indicator 
for  change  in  biodiversity.  It  does  not  reflect  individual  species 
responses  but  represents  the  average  response  of  the  total  set 
of  original  species  relative  to  their  abundance  in  undisturbed 
ecosystems  [72].  This  indicator  does  not  cover  all  aspects  of  the 
complex  concept  of  biodiversity,  but  it  can  be  used  appropriately 
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Fig.  1.  Overview  of  the  bioethanol  chain  and  the  environmental  impacts  assessed,  the  indicators  used  and  the  models  applied.  The  bottom  boxes  indicate  the  type  of  result 
that  is  generated  from  the  assessments. 


to  assess  changes  in  biodiversity  due  to  changes  in  land  use  for 
bioenergy  crops.  It  was  successfully  applied  in  several  global  en 
regional  studies  concerning  changes  in  biodiversity  [50,73,74]. 
The  MSA  values  for  several  land  cover  types  are  (included  in 
Table  A1.10  of  the  Annex)  was  linked  to  a  map  of  current  land 
use.  However,  as  this  land  use  map  does  not  account  for  diversity 
within  a  land  use  type,  management  practices  or  locally  available 
biodiversity,  the  MSA  provides  little  differentiation  for  agricultural 
land  in  the  Netherlands.  For  all  areas  outside  the  HNV  areas,  the 
MSA  indicator  is  applied  to  assess  the  impact  of  LUC  to  bioenergy 
crops  on  biodiversity. 

2.5.  Overview  methods  and  integration  of  results 

Fig.  1  shows  an  overview  of  the  bioenergy  supply  chain,  the 
impacts  that  are  assessed,  the  methods  that  are  applied  and  the 
type  of  results  that  are  produced  in  this  study. 

For  an  integral  picture,  the  GIS  maps  layers  of  the  potential 
environmental  impacts  have  been  combined.  As  all  impacts  are 
measured  in  different  units  and  scales,  combining  impacts  requires 
standardisation  of  the  GIS  map  layers.  For  all  indicators,  maximum 
standardisation  was  applied  in  which  all  negative  impacts  were 
translated  to  a  scale  from  0  to  -10  and  all  positive  effects  to  a  value 
between  0  and  +1 0.  In  this  standardisation  process,  impacts  equal  to 


the  impact  of  current  land  use  were  set  to  ‘O’,  all  increases  in  impacts 
compared  to  current  land  use  were  assumed  to  be  negative  and  all 
decreases  in  impacts  were  assumed  to  be  positive.  In  order  to  pre¬ 
vent  a  skewed  distribution  in  standardised  values  due  to  extreme 
results  of  the  indicators,  the  maxima  were  set  on  2  times  the  stan¬ 
dard  deviation  from  the  average.  In  order  to  compare  the  results 
of  Miscanthus  and  sugar  beet,  the  same  standardisation  has  been 
applied  for  the  results  of  both  crops.  In  this  integrated  assessment, 
it  is  assumed  that  all  environmental  impacts  are  equally  important. 

3.  Results 

3. 1 .  GHG  during  lifecycle 

The  GHG  emissions  over  the  production  chain  of  bioethanol  are 
substantial  (see  Fig.  2).  The  results  found  in  this  assessment  are 
comparable  with  the  figures  found  in  previous  studies  [75-78]. 
Bioethanol  production  of  Miscanthus  causes  fewer  GHG  emissions 
than  first  generation  biofuels.  As  thermal  energy  requirements  for 
2nd  generation  conversion  are  assumed  to  be  fuelled  by  part  of  the 
feedstock,  the  process  requires  less  fossil  energy.  Fig.  2  shows  the 
relative  importance  of  GHG  emissions  due  to  changes  in  SOC  com¬ 
pared  to  the  emission  during  the  lifecycle.  The  uncertainty  bars 
indicate  the  spatial  variation  in  SOC  change  for  a  20  year  period.  As 
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Fig.  2.  GHG  emissions  due  to  SOC  changes  and  during  lifecycle  of  bioethanol  production  from  Miscanthus  and  sugar  beet.  For  comparison,  the  average  GHG  emission  of 
petrol  over  the  lifecycle  is  depicted. 


the  emission  reduction  potential  of  a  bioenergy  production  chain  is 
highly  related  to  the  location  of  production  this  needs  to  be  spatially 
assessed  (see  Fig.  3). 

3.2.  GHG  emissions  due  to  LUC 

Fig.  4  shows  the  changes  in  GHG  emissions  due  to  LUC  from 
current  land  use  to  Miscanthus  and  sugar  beet.  This  includes  the 
changes  in  carbon  stocks  (Fig.  3),  the  changes  in  carbon  fluxes  from 
organic  soils  and  N20  emissions  (Fig.  A2.2  Annex). 

When  grassland  is  converted  to  Miscanthus,  a  carbon  loss  due  to 
land  clearing  (biomass  carbon,  up  to  5.4  ton  C  ha-1 )  and  soil  prepa¬ 
ration  prior  to  Miscanthus  planting  (soil  carbon,  ItonCha-1  on 
sandy  soils,  2.5  ton  C  ha-1  on  clay  soils)  is  compensated  by  relatively 
large  above  and  below  ground  biomass  stock  (up  to  1 7.5  ton  C  ha-1 ) 
that  Miscanthus  achieves  compared  to  grass.  When  cropland  is 
converted  to  Miscanthus,  relatively  large  amounts  of  carbon  are 
sequestered  both  in  biomass  and  soil  carbon  (8.5  ton  C  ha-1  on 
sandy  soils,  21.5  ton  ha-1  on  clay  soils).  For  Miscanthus,  the  GHG 
emission  reduction  through  carbon  sequestration  varies  between 
0  to  -78  kgGJethanol-1.  For  sugar  beet,  the  GHG  emissions  due  to 
changes  in  carbon  stocks  vary  between  0  and  54kgGJethanol-1. 
These  figures  are  in  line  with  the  studies  of  Foereid  et  al.  [79]  and 
Schneckenberger  and  Kuzakov  [80]  on  sequestration  potential  of 
Miscanthus  and  with  the  study  of  Kuikman  [35]  on  SOC  in  the 
Netherlands. 

LUC  related  N20  emissions  contribute  significantly  to  the  GHG 
emissions.  Conversion  of  current  land  use  to  Miscanthus  gener¬ 
ally  decreases  N20  emissions.  Especially  when  pastures  (generally 
high  fertiliser  and  manure  application  levels)  are  converted  to 
Miscanthus,  N20  emissions  decrease  significantly  (±1.7tonC02- 
eqha-1  y-1).  For  Miscanthus,  the  GHG  emission  reduction  though 
LUC  related  N20  emission  reduction  varies  between  0  and  -21  kg 
C02  eq  GJethanol-1 .  LUC  to  sugar  beet  cause  an  increase  of  N20  emis¬ 
sions  in  almost  all  areas,  with  highest  increases  on  arable  land. 
This  is  mainly  caused  by  the  higher  fertiliser  application  levels  and 


the  increased  emissions  from  crop  residues  (sugar  beet  leaves  and 
crowns).  For  sugar  beet  ethanol,  the  LUC  related  N2  O  emissions  con¬ 
tribute  0-15  kg  GJethanol-1  to  the  GHG  balance  (see  Fig.  A2.2  Annex). 

Because  LUC  could  cause  an  increase  in  one  type  of  GHG  but 
a  decrease  in  another  GHG,  changes  could  balance  each  other  out 
(for  example:  leaving  the  residues  of  sugar  beet  on  the  field  will 
increase  SOC  but  it  will  also  increase  N20  emissions  from  residues). 
In  general,  LUC  towards  Miscanthus  will  reduce  GHG  emissions, 
especially  conversion  of  arable  land.  LUC  to  sugar  beet  will  gen¬ 
erally  cause  a  net  increase  in  GHG  emissions.  The  amount  of  GHG 
emissions  is  relatively  high  in  areas  now  in  use  for  pastures.  For 
Miscanthus,  the  GHG  emission  reduction  compared  to  current  land 
use  varies  between  0  to  -159  kg  C02-eqGJethanol-1  (-15tonC02- 
eqha-1).  When  arable  land  is  converted  to  sugar  beet,  relatively 
small  changes  in  GHG  emissions  occur.  The  GHG  emissions  due  to 
LUC  to  sugar  beet  range  from  0  to  148  kgGJethanol-1  (-1  tonC02- 
eqha-1). 

3.3.  Soil  erosion 

The  risk  on  erosion  (in  ton  hay-1 )  for  current  land  use  is  shown 
in  Fig.  A2.5  in  the  Annex  and  the  changes  in  the  risk  on  erosion  due 
to  LUC  are  depicted  in  Fig.  5.  The  figures  are  based  on  the  parameters 
for  April  which  is  the  most  sensitive  period  of  the  year  for  wind 
erosion  due  to  low  crop  coverage,  decreased  soil  moisture  levels 
and  relatively  high  wind  velocities  (See  Annex  2.2). 

Although  the  overall  risk  on  erosion  is  relatively  small  (max 
<6  ton  hay-1),  this  could  still  have  severe  negative  effects.  In  the 
current  situation,  light  sandy  soils  in  the  south  east  of  the  region 
(veenkolonien)  are  the  areas  most  at  risk  for  soil  erosion.  In  the 
areas  prone  to  soil  erosion,  relative  large  changes  in  erosion  risk 
occur.  When  arable  land  is  converted  to  Miscanthus,  erosion  risk  is 
significantly  reduced.  Although  it  is  assumed  that  grass  is  renewed 
every  10  years  and  Miscanthus  every  20  years,  risk  on  erosion 
increases  when  pastures  are  converted  to  Miscanthus.  This  is 
related  to  the  planting  date  of  grass  which  is  generally  seeded  in 
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Fig.  3.  (a)  A  SOC  when  current  land  use  is  converted  to  Miscanthus.  (b)  A  SOC  when  current  land  use  is  converted  to  sugar  beet. 


September.  By  the  time  the  soil  gets  prone  to  erosion  in  spring, 
grass  already  covers  the  soil  to  a  large  extent.  Miscanthus,  however, 
needs  to  be  planted  in  spring  and  has  a  slow  crop  development  in 
the  initial  growth  stage.  Therefore,  in  the  year  of  establishment  the 
risk  on  erosion  is  relatively  high.  These  findings  are  in  line  with 
Kort  et  al.  [81].  The  risk  on  erosion  increases  up  to  3  ton  ha-1  y-1 
when  sugar  beet  replaces  crop  rotations.  When  sugar  beet  replaces 
pastures,  the  risk  on  erosion  can  even  increase  up  to  9  ton  ha-1  y_1 
on  sandy  soils. 

3.4.  Water 

3.4.1.  Water  consumption 

In  the  current  situation,  evapotranspiration  of  crop  rotation 
and  pasture  land  exceed  the  effective  precipitation  in  the  months 


April  to  September  (see  Fig.  A2.6  of  the  Annex).  Therefore,  a 
cumulative  deficiency  occurs  during  these  months.  As  the  effec¬ 
tive  precipitation  exceeds  the  evapotranspiration  considerably 
during  the  rest  of  the  year,  all  water  shortages  are  eventually 
replenished.  However,  temporary  shortages  during  the  summer 
could  cause  damages  to  agricultural  production  or  natural  areas. 
The  evapotranspiration  of  Miscanthus  exceeds  the  evapotran¬ 
spiration  of  both  pastures  and  rotations  of  arable  crops  from 
May  to  September.  During  July  and  August  the  evapotranspi¬ 
ration  of  sugar  beet  exceeds  the  average  evapotranspiration  of 
rotations  and  pastures.  Therefore,  LUC  to  bioenergy  crops  could 
cause  higher  water  consumption  rates  during  summer.  As  water 
tables  are  regulated  artificially  in  the  Netherlands,  more  water 
should  be  supplied  to  the  areas  with  higher  water  require¬ 
ments.  The  developments  in  crop  specific  evapotranspiration 


F.  van  dev  Hilst  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  16  (2012)  2053-2069 


2061 


Miscanthus 

A  GHG  (ton  C02-eq/ha) 

<-15.000 
-15000- -10000 
-1 0000  -  -500^^]  Built-Up 


]  -5000  -  0 

□  >0 


□  Nature 


Water 


0  10  20  40  Km 

I _ i _ l _ l _ I _ i _ i _ l _ I 


0  10  20  40  Km 

I _ l _ l _ l _ I _ l _ l _ i _ I 


Sugar  beet 

A  GHG  (ton  C02-eq/ha) 

<-  0.00 

0  -  5000  _ 

j  Built-up 
5000-  10000!^^ 

Nature 
Water 


10000-21 
>  20000 


000^^^ 


Fig.  4.  (a)  A  GHG  when  current  land  use  is  converted  to  Miscanthus.  (b)  A  GHG  when  current  land  use  is  converted  to  sugar  beet. 


over  time  in  relation  to  the  effective  precipitation  is  depicted  in 
Fig.  A2.6  of  the  Annex. 

The  spatial  variation  in  water  deficits  during  the  months 
April-August  caused  by  evapotranspiration  of  current  land  use  is 
depicted  in  Fig.  A2.7  of  the  Annex.  The  largest  water  shortages  occur 
in  the  western  part,  as  precipitation  is  relatively  low  in  this  area 
compared  to  the  rest  of  the  region.  The  change  in  water  shortages 
when  current  land  use  is  converted  to  Miscanthus,  are  depicted  in 
Fig.  6a.  The  spatial  variation  in  change  in  water  deficits  due  to  the 
introduction  of  sugar  beet  is  depicted  in  Fig.  6b. 

In  addition  to  the  spatial  water  balance,  the  regional  WUE  was 
measured.  In  Fig.  A2.8  of  the  Annex,  the  WUE  for  Miscanthus  and 
sugar  beet  is  depicted.  Most  studies  indicate  that  the  WUE  of  C4 
crops  is  generally  higher  than  the  WUE  of  C3  crops  [48,50].  However 


in  this  study,  it  was  found  that  the  WUE  of  Miscanthus  (2.23  g  kg-1 
water)  is  lower  that  for  sugar  beet  (3.87  g  kg-1  water).  This  can 
be  explained  by  the  fact  that  the  competitive  advantage  of  the  C4 
photosynthesis  pathway  decreases  significantly  further  away  from 
the  equator.  In  addition,  as  climatic  conditions  and  agricultural 
management  are  optimal  for  sugar  beet  production,  it  achieves  rel¬ 
atively  high  WUE.  The  WUE  found  for  Miscanthus  is  in  line  with  the 
2.1  g  kg-1  found  by  Clifton  Brown  [82]  and  the  range  of  1  -9.5  g  kg-1 
provided  by  Berndes  [48].  The  regional  WUE  of  sugar  beet  is  in  line 
with  the  range  of  2.55-3.83  g  kg-1  given  by  [50]. 

3.4.2.  Water  quality 

The  Miterra  model  was  used  to  calculate  the  nitrogen  and  phos¬ 
phorous  balance  in  the  soil.  Based  on  the  N  balance  and  the  water 
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Fig.  5.  (a)  A  Wind  erosion  risk  in  kg  soil  ha-1  y-1  when  current  land  use  is  converted  to  Miscanthus.  (b)  A  Wind  erosion  risk  in  kg  soil  ha-1  y-1  when  current  land  use  is 
converted  to  sugar  beet. 


surplus,  the  N03_  concentration  in  the  ground  water  was  calcu¬ 
lated.  The  maps  of  the  N  and  P  balance  are  included  in  Annex 
2.3.  As  N  and  P  fertiliser  application  levels  are  closely  linked,  the 
spatial  patterns  of  changes  in  N  and  P  balances  are  quite  similar. 
Because  of  the  large  differences  between  the  current  fertiliser  and 
manure  input  on  pasture  and  the  fertiliser  requirements  of  Mis¬ 
canthus,  the  N  and  P  surplus  will  be  significantly  reduced  when 
pasture  area  is  converted  to  Miscanthus.  Also  in  the  areas  where 
sugar  beet  substitutes  grassland  N  and  P  surpluses  are  reduced. 
However,  when  sugar  beet  substitutes  a  rotation  of  crops,  N  and  P 
surpluses  increase.  The  high  increases  of  N  and  P  occur  when  rota¬ 
tions  with  high  proportions  of  cereals  are  converted  to  sugar  beet. 
This  is  mainly  caused  by  the  relatively  high  fertiliser  application 
rates  of  sugar  beet  compared  to  other  rotation  crops. 


The  differences  in  patterns  between  the  N  balance  and  N03_ 
concentrations  can  be  explained  by  the  differences  in  leaching  rates 
(lower  for  grassland,  peat  soils  and  lower  precipitation  surpluses). 
The  N03-  concentration  in  mg  N03  l-1  is  shown  in  Fig.  7.  The  cur¬ 
rent  N03-  levels  are  72  mgl-1  on  average  and  75%  of  the  zip  code 
districts  exceed  the  standards  of  50gl_1  of  the  Nitrate  directive 
[83].  In  general,  the  N03-  concentration  decreases  when  current 
land  use  is  converted  to  Miscanthus.  The  average  change  in  N03- 
concentrations  is  -53  mgH1  and  varies  between  0  in  areas  with  a 
large  fraction  of  arable  land  and  - 1 87  mg  1_1  in  zip  code  areas  with 
a  high  proportion  of  grassland.  When  pasture  area  is  converted  to 
sugar  beet,  a  decrease  in  N03-  concentration  (of  max  -75mgl_1) 
occurs.  When  arable  land  is  converted  to  sugar  beet,  an  average 
increase  in  N03-  concentrations  of  16.5  mg l-1  is  simulated. 
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Fig.  6.  (a)  A  Water  depletion  during  summer  (mm)  when  current  land  use  is  converted  to  Miscanthus.  (b)  A  Water  depletion  during  summer  (mm)  when  current  land  use 
is  converted  to  sugar  beet. 


3.5.  Biodiversity 

In  Fig.  8,  the  HNV  areas  were  biodiversity  is  at  risk  when  current 
land  use  is  converted  to  bioenergy  crops  are  mapped.  It  shows  that 
in  almost  the  entire  pasture  area  in  the  north  of  the  Netherlands 
which  comprises  high  densities  of  important  species  (endangered 
bird  species  and  large  shares  of  European  populations  of  meadow 
and/or  wintering  birds),  biodiversity  is  at  risk  when  current  land 
use  is  converted  to  either  Miscanthus  or  sugar  beet. 

The  HNV  areas  where  biodiversity  is  at  risk  when  converted  to 
sugar  beet  (532  x  103  ha,  of  which  193  x  103  ha  at  ‘high  risk’)  are 
more  extensive  than  for  conversion  to  Miscanthus  (473  103  ha,  of 
which  91  x  103  ha  at  ‘high  risk’).  The  areas  where  LUC  to  bioenergy 
crops  is  expected  to  cause  a  considerable  risk  of  biodiversity 
loss,  are  considered  to  be  areas  where  energy  crops  should  be 


introduced  with  care.  The  areas  with  a  high  risk  on  biodiversity 
loss  are  considered  to  be  no-go  areas  for  bioenergy  crops.  In  some 
areas  (59  x  103  ha)  the  introduction  of  Miscanthus  could  have  both 
a  positive  and  a  negative  effect  on  biodiversity.  These  are  mainly 
the  heterogeneous  areas  and  the  extensive  managed  arable  areas 
in  the  east  of  the  region. 

For  the  areas  not  included  in  the  HNV  areas,  the  change  in  bio¬ 
diversity  is  expressed  in  the  change  in  MSA  value.  This  mapped  in 
Fig.  Al .  1 1  a  and  b  of  the  Annex.  The  change  in  MSA  value  indicates  an 
increase  in  biodiversity  when  pasture  land  is  converted  to  Miscant¬ 
hus  and  especially  when  arable  land  is  converted  to  Miscanthus. 
The  MSA  values  indicate  a  decrease  in  biodiversity  when  pasture 
land  is  converted  to  sugar  beet  cultivation  and  a  status  quo  in  bio¬ 
diversity  when  sugar  beet  is  cultivated  in  arable  areas.  The  results 
show  that  sugar  beet  will  induce  a  risk  of  biodiversity  loss  both 
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Fig.  7.  (a)  A  N03_  concentration  when  current  land  is  converted  to  Miscanthus.  (b)  A  N03_  concentration  when  current  land  is  converted  to  sugar  beet. 


within  and  outside  the  HNV  areas  and  that  the  only  areas  where  no 
negative  effects  on  biodiversity  occur  are  the  areas  already  in  use 
for  intensively  managed  arable  crops. 


3.6.  Integrated  results 

In  Fig.  9a  and  b  the  integrated  results  of  all  environmental 
impacts  are  spatially  depicted  for  both  Miscanthus  and  sugar  beet. 
The  results  are  expressed  as  an  average  score  on  a  standardised 
scale  from  -10  to  +10.  In  general,  Miscanthus  has  better  overall 
scores  (mean  =  2.33)  than  sugar  beet  (mean  = -0.85).  The  results 
discussed  in  previous  sections  show  a  wide  distribution  of  scores 
for  the  individual  impacts.  As  the  integrated  results  show  average 
scores  close  to  0,  it  illustrates  that  there  are  several  trade-offs 


between  impacts  and  that  positive  and  negative  scores  balance 
each  other  out  to  some  extent. 

There  are  no  areas  where  only  positive  impacts  occur  when 
current  land  use  is  changed  to  sugar  beet  or  Miscanthus.  For  both 
crops,  the  western  wet  pasture  area  and  the  peaty  south  east  area 
of  Drenthe  appear  to  be  the  areas  with  most  negative  impacts.  The 
conversion  of  pasture  land  to  arable  crops  on  (reclaimed)  peat  areas 
has  negative  impact  on  carbon  stocks,  risk  on  erosion  and  biodiver¬ 
sity.  These  negative  impacts  are  relatively  small  for  Miscanthus  but 
quite  severe  for  sugar  beet.  The  impacts  in  these  areas  would  even 
be  more  severe  in  case  the  pasture  areas  were  not  as  intensively 
managed  as  they  are  today.  Sugar  beet  has  also  a  negative  envi¬ 
ronmental  performance  in  the  eastern  part  of  the  region.  This  area 
is  characterised  by  arable  land  on  sandy  soils.  As  sugar  beet  is  a 
more  intensively  managed  crop  than  average  arable  crop  rotations, 
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Fig.  8.  (a)  Risk  of  biodiversity  loss  when  current  land  use  is  converted  to  Miscanthus  based  on  HNV  cluster  indicator,  (b)  Risk  of  biodiversity  loss  when  current  land  use  is 
converted  to  sugar  beet  based  on  HNV  cluster  indicator. 


negative  impacts  on  the  N03_  concentration,  P  balance,  soil  erosion 
and  also  biodiversity  are  relatively  large.  For  the  same  reasons, 
Miscanthus  (extensively  managed  crop)  scores  relatively  well  in 
this  area.  Sugar  beet  has  the  least  negative  impacts  in  the  northern 
pasture  areas  on  clay  soils  and  negative  impacts  are  partly  balanced 
out  by  reduced  risk  of  nutrient  leakage  and  water  consumption. 

Comparing  the  integrated  results  of  the  environmental  impacts 
assessed  in  this  study  with  the  results  of  the  economic  performance 
of  sugar  beet  and  Miscanthus  in  this  region  assessed  in  the  study  of 
Van  der  Hilst  et  al.  [22],  show  that  the  areas  with  best  environmen¬ 
tal  performance  are  also  the  areas  with  the  least  cost  of  production 
(in  €GJbiomass-1 ).  These  are  the  areas  currently  in  use  for  intensive 


agriculture.  However,  in  most  of  these  areas  bioenergy  crops  can¬ 
not  compete  with  current  land  use  as  conventional  crops  have  a 
better  economic  performance  on  these  suitable  soils.  Based  on  the 
results  of  the  environmental  impacts  and  economic  performance, 
an  area  of  24  kha  can  be  identified  as  most  promising  for  Miscant¬ 
hus  as  in  these  areas  it  could  compete  with  current  land  use,  it  has 
relatively  low  production  costs  and  has  the  least  negative  and  also 
positive  environmental  impacts.  This  well  exceeds  the  proportion 
of  agricultural  land  that  could  become  available  for  bioenergy  crop 
production  in  201 5  (9.6  kha  arable  land  and  1 .5  kha  pasture  land  in 
the  considered  region)  as  indicated  in  the  Refuel  study  [26].  As  the 
study  of  Hilst  et  al.  [22]  showed  that  sugar  beet  cannot  compete 
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Fig.  9.  (a)  Integrated  result  of  environmental  impact  of  a  shift  from  current  land  use  to  Miscanthus  based  on  standardised  maps  of  individual  environmental  impacts,  (b) 
Integrated  result  of  environmental  impact  of  a  shift  from  current  land  use  to  sugar  beet  based  on  standardised  maps  of  individual  environmental  impacts. 


with  current  land  use  and  has  relative  high  production  costs  in  the 
whole  agricultural  area  in  the  North  of  the  Netherlands  and  this 
study  shows  mainly  negative  environmental  impacts,  it  is  assumed 
that  sugar  beet  is  not  to  be  a  promising  bioenergy  crop  in  this  region. 

4.  Discussion  and  conclusions 

Preferably,  impacts  would  be  assessed  on  most  detailed  spa¬ 
tial  level.  However,  several  input  parameters  of  the  Miterra  model 
are  only  available  on  a  postal  code  level.  In  addition,  crop  areas 
on  a  zip  code  level  provides  a  proper  indication  of  the  crop  mix 
in  rotation,  whereas  a  land  cover  recording  of  a  single  year  on 
a  detailed  grid  level  ignores  the  effect  of  rotating  crops.  Further¬ 
more,  as  there  are  uncertainties  in  several  key  input  parameters 


such  as  crop  management,  leakage  fraction  per  soil  type,  water 
surplus  etc.,  impact  assessment  on  a  detailed  grid  level  could 
give  a  too  optimistic  impression  on  the  accuracy  of  the  results. 
There  are  some  drawbacks  of  the  Miterra  model:  The  fertiliser- 
derived  N20  emissions  are  calculated  in  accordance  with  the  Tier 
1  approach  of  the  IPCC.  However,  as  indicated  by  Smeets  et  al.  [84] 
and  Leschen  et  al.  [85]  this  neglects  the  variability  of  emissions 
due  to  differences  in  environmental  conditions,  management  and 
crop  type.  The  Miterra  Model  includes  only  regional  averages  of 
climate  data  and  in  the  calculation  of  leaching  of  nutrients  and 
decomposition  of  organic  matter,  water  surplus  of  precipitation 
was  included  but  water  table  levels  were  neglected.  Despite  these 
shortcomings,  the  Miterra  model  was  able  to  combine  several  spa¬ 
tial  datasets  in  a  consistent  way  and  therefore  provides  excellent 
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information  on  the  relative  differences  in  impacts  on  a  spatially 
detailed  level. 

As  indicated  by  Hoefnagels  et  al.  [76],  Cherubini  et  al.  [88]  and 
van  Dam  et  al.  [89],  GHG  balances  of  biofuels  are  very  sensitive 
for  assumptions  regarding  system  boundaries,  allocation  methods 
and  input  parameters.  Therefore,  the  results  of  the  GHG  assessment 
should  be  interpreted  with  care. 

High  soil  losses  due  to  wind  erosion  occur  only  incidentally 
when  all  relevant  conditions  are  disadvantageous  at  the  same  time. 
By  using  an  average  monthly  climate  factor  in  the  WEQ  no  accu¬ 
rate  indication  is  provided  on  how  often  erosion  thresholds  are 
met.  Although  the  WEQ.  is  not  the  most  sophisticated  approach  to 
assess  the  risk  on  erosion,  it  does  provide  a  method  to  discriminate 
between  crop  covers  and  to  assess  the  spatial  variability  in  the  ero¬ 
sion  risk  on  a  regional  level.  It  is  recommended  to  use  actual  field 
measurements  and  more  exact  models  to  spatially  assess  erosion 
risk  on  a  regional  level  in  order  to  achieve  more  accurate  figures  on 
the  exact  amount  of  soil  loss. 

The  cumulative  water  deficit  is  assessed  taking  the 
spatial-temporal  precipitation  and  evapotranspiration  levels 
into  account.  It  neglects  however  groundwater  level,  flow  schemes 
of  the  water  basin  and  the  extensive  water  table  regulation 
system.  Therefore,  the  spatial  water  balance  provides  too  little 
information  to  assess  where  actual  droughts  or  water  surpluses 
will  occur.  In  order  to  assess  the  actual  effect  on  water  tables 
it  is  recommended  that  more  advanced  hydrologic  models  are 
applied. 

In  this  study  the  HNV  and  MSA  indicators  were  applied  to  iden¬ 
tify  areas  with  a  high  risk  on  biodiversity  loss.  The  HNV  approach 
enables  to  exclude  the  most  vulnerable  areas  from  a  biodiversity 
perspective,  which  in  line  with  the  bioenergy  sustainability  cri¬ 
teria  such  as  developed  by  RSB,  RSPO,  NTA8080.  However,  due 
to  a  lack  of  information  on  the  response  of  species  on  the  intro¬ 
duction  of  bioenergy  crops,  no  quantitative  assessment  could  be 
made.  The  MSA  indicator  on  the  other  hand,  enables  to  quantify 
the  overall  effect  of  a  shift  to  bioenergy  crops  on  biodiversity  but 
because  the  MSA  values  does  not  account  for  diversity  within  a 
land  use  type,  management  practices  or  locally  available  biodiver¬ 
sity,  the  MSA  provides  little  differentiation  for  agricultural  land  in 
the  Netherlands. 

In  order  to  assess  the  impact  on  biodiversity  in  a  quantitative 
way,  we  conclude  in  line  with  Rowe  [18]  that  more  empirical  data 
are  required  on  the  impact  of  specific  bioenergy  crops  on  the  occur¬ 
rence  of  individual  species. 

The  calculations  of  the  environmental  impacts  are  all  based  on 
a  very  broad  range  of  input  parameters.  As  the  results  are  affected 
by  the  assumptions  made  on  climate,  soil,  agronomy  and  manage¬ 
ment  applied,  they  should  be  interpreted  with  care.  However,  as 
the  assumptions  are  consistently  applied,  it  provides  a  good  indi¬ 
cation  the  relative  spatial  differences  in  environmental  impacts.  As 
all  impacts  are  related  to  ecosystem  functionality,  they  are  heav¬ 
ily  interlinked  with  each  other.  In  some  cases,  impacts  reinforce 
each  other  and  in  other  cases,  trade  offs  between  impacts  occur. 
Although  it  is  hard  to  prove  causality  and  to  quantify  relations, 
it  would  be  interesting  to  assess  the  spatial  correlation  between 
several  impacts. 

All  environmental  impacts  of  Miscanthus  and  sugar  beet  cul¬ 
tivation  have  been  integrated  in  a  standardised  map  of  overall 
environmental  performance.  This  required  an  interpretation  of  the 
dose-response  relationship:  to  what  extend  is  the  caused  effect 
a  negative  or  positive  impact.  Both  the  dose-response  relation¬ 
ship  and  the  thresholds  for  actual  damage  are  expected  to  be 
spatially  heterogeneous,  time  related,  and,  in  some  cases,  con- 
tradictive.  In  addition,  cause-effect  relationships  could  be  linear, 
quadratic,  exponential,  inversely  proportional,  etc.  As  there  is 
too  little  knowledge  available  of  the  dose-response  relationship, 


a  simple  approach  of  maximum  linear  standardisation  with  no 
threshold  was  applied  in  this  study.  More  expertise  is  required 
on  the  translation  of  a  change  in  value  of  an  indicator  to  an 
actual  impact.  This  requires  not  only  knowledge  on  physical  rela¬ 
tionships  but  also  on  the  focal  points  of  policy,  the  willingness 
to  accept  and  the  different  points  of  view  of  the  stakeholders 
involved.  In  the  integration  of  environmental  impacts  all  impacts 
were  considered  to  be  equally  important,  whereas  ideally  multiple 
perspectives  from  stakeholder  on  the  hierarchy  of  environmen¬ 
tal  impacts  would  be  included  in  a  spatially  explicit  Multi  Criteria 
Analysis  (MCA). 

This  first  explorative  step  could  be  further  developed  to 
an  integrated  assessment  method  by  means  of  a  stakeholder 
assessment  and  more  advanced  standardisation  and  weighting 
methods. 

The  integrated  assessment  shows  that  there  are  no  areas  in 
the  North  of  the  Netherlands  where  only  positive  effects  occur 
when  bioenergy  crops  are  introduced  and  that  there  are  sev¬ 
eral  trade-offs  between  the  impacts.  In  general,  sugar  beet  causes 
relatively  many  negative  environmental  impacts  especially  in  pas¬ 
ture  areas.  In  these  areas,  the  LUC  related  GHG  emissions  are 
up  to  148  kgGJethanol_1  and  the  risk  on  soil  erosion  increases  to 
9  ton  ha-1.  Also,  there  is  a  high  risk  on  biodiversity  loss  in  these 
areas.  On  the  other  hand,  a  shift  towards  sugar  beet  could  cause 
a  decrease  of  75  mg  l-1  in  the  N03_  concentration  of  groundwa¬ 
ter  and  a  decrease  in  the  seasonal  water  deficits  with  100  mm. 
LUC  to  Miscanthus  has  less  negative  environmental  impacts  and 
in  the  areas  currently  in  use  as  arable  land  many  positive  envi¬ 
ronmental  impacts  are  expected  to  occur.  In  these  areas,  the  LUC 
related  GHG  emissions  are  reduced  by  -159  kg  GJethanol_1,  the 
risk  on  soil  erosion  could  be  reduced  with  4  ton  ha-1,  the  N03_ 
concentration  is  reduced  with  53mgl-1,  and  the  shift  towards 
Miscanthus  could  have  a  positive  contribution  to  biodiversity.  On 
the  other  hand,  at  some  locations  the  seasonal  water  depletion 
increases  with  150  mm  arable  land  is  converted  to  Miscanthus. 
For  both  crops,  the  western  wet  pasture  areas  appear  to  be 
the  area  with  most  negative  impacts.  It  can  be  concluded  that 
the  spatial  correlation  between  impacts  and  current  land  use  is 
higher  than  between  impacts  and  soil  type.  The  combination  of 
the  spatial  distribution  of  the  environmental  impacts  and  the 
economic  performance  of  sugar  beet  and  Miscanthus  assessed 
in  the  study  of  Van  der  Hilst  et  al.  [22]  shows  that  the  areas 
with  best  environmental  performance  are  also  the  areas  with  the 
least  cost  of  production  (in  €GJbiomass_1 ).  Considering  the  large 
amount  of  negative  impacts  on  areas  currently  in  use  as  pastures, 
it  could  be  recommended  to  exclude  these  areas  for  bioenergy 
crop  production.  In  addition,  as  sugar  beet  causes  relative  many 
negative  impacts,  large  scale  introduction  of  this  crop  for  bioen¬ 
ergy  purposes  should  be  carefully  considered  and  implemented 
wisely. 

This  assessment  has  resulted  in  understanding  in  which  areas 
the  introduction  of  bioenergy  is  likely  to  cause  severe  environmen¬ 
tal  impacts  and  in  which  areas  little  negative  and  even  positive 
effects  are  caused.  This  study  well  illustrates  the  high  spatial  vari¬ 
ability  of  environmental  impacts  of  potential  bioenergy  chains  and 
therefore  it  stresses  the  importance  to  assess  the  sustainability  of 
bioenergy  supply  chains  in  an  integrated  and  spatially  explicit  way. 
This  study  provides  a  set  of  methodologies  to  assess  the  spatial  vari¬ 
ation  of  environmental  impacts  of  bioenergy  chains.  Therefore,  it 
contributes  to  identify  promising  locations  for  bioenergy  produc¬ 
tion  from  an  environmental  point  of  view.  Because  of  the  generic 
characteristics,  the  methodologies  and  tools  applied  for  this  study 
can  also  be  used  for  other  case  studies.  It  should  be  noted  that 
this  type  of  spatial  explicit  regional  assessment  on  environmental 
impact  of  bioenergy  chains  cannot  replace  site  specific  environ¬ 
mental  impact  assessment. 
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